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ABSTRACT: Here, we studied the genetic diversity of native fowls in Laos by analyzing a mitochondrial DNA (mtDNA) sequence 
polymorphism. A 546-bp fragment of the mtDNA D-loop region was sequenced in 129 chickens from the areas of Vientiane, Luang 
Prabang and Pakse. In total, 29 haplotypes were identified and formed five clades. Haplotype diversity and nucleotide diversity of the 
native fowls in Laos were 0.85536±0.0172 and 0.010158±0. 005555, respectively. Although the Laotian native fowls were distributed 
across five clades, most of them were clustered in two main clades (A and B), which were originated in China. The other haplotypes 
were contained in clades D, F, and I, which originated from continental southeast Asia. These results suggest that multiple maternal 
lineages were involved in the origin of domestic chicken in Laos. Moreover, there appear to be at least two maternal lineages, one from 
China and the other from the southeast Asian continent. (Key Words: Native Fowl, Genetic Variability, mtDNA, Maternal Origin) 



INTRODUCTION 

Laos is a tropical landlocked country that shares borders 
with Thailand, Cambodia, Vietnam, China, and Myanmar, 
and is located between latitudes 14 and 23 and longitudes 
100 and 108. The country's borders with Thailand and 
Myanmar are formed by the Mekong River, whose source is 
in Quinghai Province, China. The terrain in over 70% of the 
country is composed of mountains and plateaus. There are 
many kinds of native fowls in Asia. Most of them have not 
been improved and have lower productive performance than 
the improved western breeds. In Laos, there are many kinds 
of native fowls, and farmers rear small flocks of native 
chickens as an important protein source. There are also 
many improved breeds that have been introduced from 
foreign countries, mainly Thailand. Pure native fowl in 
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Laos are gradually decreasing. Because the native chicken 
populations possess genes that modern breeds have lost, it 
is very important to retain these genes for the future. Thus, 
we need to survey the genetic characteristics of native fowls 
and evaluate their genetic resources. However, there has 
been little research on Laotian native fowls. One of the few 
studies was done by Okamoto et al. (1996), who reported 
the gene constitution of Laotian native fowl by analyzing 
blood protein polymorphisms. 

Many researchers have reported mitochondrial DNA 
(mtDNA) polymorphisms in various animal populations to 
clarify the ancestral lineage or to compare populations. A 
number of studies have been caiTied out to investigate the 
ancestors of modern chickens by examing mtDNA 
polymorphisms. Akishinonomiya et al. (1996) suggested a 
monophyletic origin of domestic chickens from one of the 
red jungle fowl subspecies Gallus gallus gallus, and 
proposed that a single domestication event occurred in 
Thailand and adjacent countries. In contrast, Kanginakudru 
et al. (2008) reported that two other subspecies, G. g. 
spadiceus and G. g. murghi, also contributed to the lineage 
that gave rise to domestic chickens. Other repoits described 
multiple and independent domestication events in south 
China, southeast Asia and the Indian subcontinent (Liu et al.. 
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2006; Oka et al., 2007). Liu et al. (2006) revealed nine 
divergent clades related to the geographical distribution of a 
wide range of domestic chickens in Eurasian regions. Oka 
et al. (2007) also identified seven clades in Japanese native 
chickens, of which four clades are identical to four clades 
described by Liu et al. (2006). Muchadeyi et al. (2008) 
arranged these clade classification approaches by Liu et al. 
(2006) and Oka et al. (2007) to apply Zimbabwean native 
chickens analysis, and suggested that a third maternal 
lineage excluded Zimbabwean and other African chickens 
and clustered with haplotypes presumably originating from 
China. Cue et al. (2011) also apphed this strategy to 
Vietnamese native breeds and suggested that they have 
originated from multiple maternal lineages from several 
regions of China and surrounding regions. 

The present study was conducted to clarify the genetic 
diversity of Laotian native fowl populations. We also tried 
to determine the degree of shared maternal mtDNA 
haplotypes among populations of Laotian native fowls and, 
hence, reveal maternal hneages of origin. 

MATERIALS AND METHODS 

Birds 

A total of 129 native chicken samples were collected for 
this study. The native fowls were collected 54 birds from 
Vientiane, 42 from Luang Prabang, and 33 from Pakse. 

mtDNA amplification and sequencing 

Total DNAs were extracted from blood samples by 
using a standard phenol-chloroform extraction. A fragment 
of 546-bp from the mtDNA D-loop region was amplified 
using PCR. The primers used were L16750 (5'- 
AGGACTACGGCTTGAAAAGC-3'; Akishinonomiya et al., 
1994) and H522 (5'- ATGTGCCTGACCGAGGAACCAG - 
3'; Liu et al., 2006). The niunbers in the primer names 



indicate the homologous positions of the 3' end of the 
primers on the mtDNA sequence described by Desjardins 
and Morals (1990). The PCR reaction was performed using 
the GeneAmp PCR System 9700 (Applied Biosystems, CA, 
USA) and the following mixture consisted of 100 ng of 
template DNA, IxPCR reaction buffer, 4 pmol of each 
primer, 400 jimol of each dNTPs, and 1 U of exTaq 
polymerease (TaKaRa, Otsu, Japan). The thermal profile 
included an initial denaturation at 94 °C for 1 min followed 
by 30 cycles, each of which included denaturation at 94°C 
for 1 min, annealing at 60°C for 1 min, and extension at 
72°C for 1 min, with a final extension step at 72°C for 7 
min. The PCR products were isolated from 1 % agarose gels 
and purified with Mag Extractor (TOYOBO, Osaka, Japan). 
The BigDye terminator cycle sequencing kit v. 3.1 (Applied 
Biosystems) and ABl Prism 3130x1 genetic analyzer 
sequencer (Apphed Biosystems) were used for sequence 
analysis. 

Data analysis 

Multiple alignment analysis was carried out using the 
Clustal X version 2.1 computer software (Larkin et al., 
2007). The position and number of polymorphic sites and of 
corresponding haplotypes were calculated using MEGA 
v.5.2.1 (Tamura et al., 2011). Nucleotide diversity (Nei and 
Li, 1979) and haplotype diversity (Nei, 1987) were 
estimated using ARLEQUIN v.3.5.1.3 (Excoffier et al., 
2010). A median joining network (Bandelt et al., 1999) was 
constructed using NETWORK 4.611 software (Fluxus 
Technology Ltd.) to classify the haplotypes under the 
settings described by Cue et al. (2011) into nine clades, 
following Liu et al. (2006), and three clades, following Oka 
et al. (2007). The list of haplotypes used and the 
corresponding GenBank accession numbers are provided in 
Table 1. 



Table 1. Haplotype names and accession numbers of chicken mtDNA sequences used in this study 



Haplotype GenBank accession No. Reference 



Al -AlO 




This study 


Bl -B12 




This study 


D1-D4 




This study 


F1-F2 




This study 


I 




This study 


Liu_Al 


AB 11 4069 


Liu et al. (2006) haplotype Al 


Liu_Bl 


AB007744 


Liu et al. (2006) haplotype Bl 


Liu_Cl 


AB 114070 


Liu et al. (2006) haplotype CI 


Liu_Dl 


AY588636 


Liu et al. (2006) haplotype Dl 


Liu_El 


AB 114076 


Liu et al. (2006) haplotype El 


Liu_Fl 


AF5 12285 


Liu et al. (2006) haplotype Fl 


Liu_Gl 


AF512288 


Liu et al. (2006) haplotype Gl 


Liu_Hl 


D82904 


Liu et al. (2006) haplotype HI 


Liu_Il 


AB009434 


Liu et al. (2006) haplotype \l 


Oka_D6 


AB268535 


Oka et al. (2007) haplotype D6 


Oka_Gl 


AB268545 


Oka et al. (2007) haplotype Gl 


Oka_Fl 


AB268543 


Oka et al. (2007) haplotype Fl 
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Table 2. The identified liaplotypes with the mitochondrial D-loop sequcnec polymorphisms in Laotian native fowl populations 

Nucleotide substitution 



1 1 1 1 12222222222222222222222222222233333333333333334444 
679991111222223333334444455668899900111124455566791444 
770890279057890346790236946151356926035722724537017679 



Al 


CATCTCGTCTTCCCTCCTAACGTCATTCCAATCTCTCCCATAACTTCTTCCCTT 


A2 




r 












A3 








G T 








A4 










T 






A5 




c 










G 


A6 














G 


A7 






A 










A8 














T 


A9 






c 










AlO 














c 


Bl 


T 




A 


c 


T 


T 




B2 


T 




A 


c 


T T 


T 




B3 


T 




A 


T C 


T 


T 




B4 


T 




A 


c 


C T 


T 




B5 


T 




A 


C T 


T 


T 




B6 


T 




. . A . 


. . C 


. . T 


. . . T T 




B7 


T 




. . A . 


. . C 


. . T 


. . . . T 


G 


B8 


T 




, . . A . 


. . C G 


. . T 


. . . . T 




B9 


T 




, . . A . 


. . C 


..T CT. 


. . . . T 




BIO 


T 




. . . A . 




. . T T . 


. . . . T 


G 


Bll 


T 




, . . A . 


. . C 


. . T T . 


. . . . T 


G 


B12 


T 




, . . A . 


. . C 


.CT C... 


. . . . T 




Dl 


T 




, . T . . 


. . C 


.C...CT.G.... 


. CT . . . . G . 




D2 


T 






. . C 


.C...CT.G.... 


. C T 


. . . . C 


D3 


T 






. C C 


.C...CT.G.... 


. C T 




D4 


T 




, . . A . 


. . C 


.C...CT.G.... 


. C T 




Fl 


T 






..C TC... 


.C..CCT 


. . T . T C . . . 


T . . 


F2 


T 






..C.T...TC... 


.CT.CCT 


. . T . T C . . . 


T . . 


I 


T . 




, . . A . 


..C.T G.. 


.CT G.... 


. . T 


. . C T . . T . . . . 



Dot (.) represents the identical nucleotide with the haplotype Al. 



RESULTS AND DISCUSSION 

Genetic diversity and haplotype classification 

A total number of 37 variable sites were identified 
(Table 2). No insertions or deletions were detected in our 
sequences. The calculated number of polymorphic sites and 
haplotypic diversity of three populations of Laotian natives 
are shown in Table 3 . The haplotypic diversity ranged from 
0.5183 (Luang Prabang) to 0.8798 (Vientiane), although the 
differences between populations were small. Laotian native 
fowls showed a higher value of haplotypic diversity than 
did Zimbabwean (0.29 to 0.78: Muchadeyi et al., 2008), 



Vietnamese (0.62 to 0.94: Cue et al., 2011) or Korean (0.59 
to 0.82: Hoque et al., 2013) fowls. The nucleotide diversity 
observed in this study ranged from 0.0007706 (Luang 
Prabang) to 0.014519 (Pakse) and was almost the same 
value as that of Korean native chickens (Hoque et al., 2013). 
These results suggest that Laotian native chicken 
populations have relatively high genetic diversity. 

Network analysis and haplotype distribution 

The 29 haplotypes observed in the Laotian native fowls 
were clustered into five clades (Figure 1). The haplotypes 
Al, Bl, Dl, Fl, and I in this study were the same as the 



Table 3. Polymorphic sites, haplotype and nucleotide diversity of three Laotian native fowl populations 



Locality 



No. of 
polymorphic sites 



No. of 
haplotypes 



Haplot5'pe diversity (+SD) 



Nucleotide diversity (+SD) 



Vientiane 54 

Luang Prabang 42 

Pakse 33 

Total 129 



23 
10 
26 
37 



20 
7 
10 
29 



0.8798±0.0268 
0.8153±0.0296 
0.8466±0.0350 
0.8536+0.0172 



0.008957+0.005033 
0.007706±0.004448 
0.014519+0.007823 
0.010158+0.005555 
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Clade A 



A02 A03 .AOS ^» 

^ ^. P i B12 



Clade B 



Clade F 

F01 


UuE1( 


UuFI Cl 






F02« 


UuGI 





• -AID 

A04 



Clade D 





Clade I 



• Vientiane Luang Prabang • Pakse 

• median vectors used in connecting indirectly related haplotypes 

Figure 1. Median networli profile of tfie mtDNA Dloop liaplotypes observed in the present study. Data are merged with sequences of 
major haplotypes described from Liu et al. (2006) and Oka et al. (2007). The circle size con'esponds to haplotype frequency. 

partial sequence of each haplotype Al, Bl, Dl, Fl, and II 
from the clades described by Liu et al. (2006). We did not 
detect the five haplotypes Cl, Dl, El, Fl, and Gl described 
by Liu et al. (2006) or the three haplotypes D6, Fl, and Gl 
described by Oka et al. (2007) in the study populations. 
Most individuals were classed in clades A and B whereas 
clades D, F, and I consisted of only a small number of 
individuals (Table 4). 

The formation of the network profile of Laotian native 
fowl populations was similar to that of Vietnamese native 
chicken breeds reported by Cue et al. (2011). However, 
clades C, E, and G which are major haplotypes in 
Vietnamese chickens, were not observed in the Laotian 
populations. Within clades A and B, the major haplotypes 
were Al and Bl, composing 59% and 46% of each clade, 
respectively. Using the skeleton of supposed regions of 
domestication, this finding suggests that two maternal 
lineages dominate in the Laotian native fowls, which 
presumably originated from Yunnan and the surrounding 
area in China (Liu et al., 2006). A total of 80% of the Pakse 
population belonged to clades D, F, and L Liu et al. (2006) 
and Oka et al. (2007) suggested that clade D has its roots in 
southeast, south and southwest China and/or the 
surrounding area (i.e., Vietnam, Thailand, Myanmar and 
India). A small number of chickens in the Pakse population 
were also distributed in clade I, originating from continental 
southeast Asia (Liu et al., 2006). Both Vientiane and Luang 
Prabang are northern areas in Laos, and they are close to 
China. However, Pakse area is located in the southern part 
of Laos is close to Vietnam. Therefore, our results are 
consistent with the geographical relationship of these 
populations. 



Table 4. 


Distribution of mtDNA D-loop haplotypes 


in three 


populations of Laotian native fowls 






Haplotype 


Vientiane 


Luang 
Prabang 


Pakse 


Total 


Al 


14 


8 


8 


30 


A2 


3 


5 


1 


9 


A3 




4 




4 


A4 




2 




2 


A5 


1 








A6 


1 








A7 




1 






A8 


1 








A9 


1 








AlO 


1 








Bl 


9 


13 


9 


31 


B2 


9 


9 


5 


23 


B3 


4 








B4 










B5 










B6 










B7 










B8 










B9 






1 




BIO 










Bll 










B12 










Dl 






4 




D2 










D3 










D4 






1 




Fl 






1 




F2 






1 




I 






2 


2 


Total 


54 


42 


33 


129 
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